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A new primary mass independent method of energy measurement has been developed by exploit- 
ing: a) the joint analysis of the shower size, obtained by EAS-TOP, and of the EAS atmospheric 
Cherenkov light lateral distribution, obtained by the QUEST array, and b) simulations based on the 
CORSIKA/QGSJET code. The method is based on the correlation between the size/energy ratio and 
the steepness of Cherenkov light lateral distribution and has been compared with a "classical" one 
based on the Cherenkov light flux at a fixed distance (175 m) from the EAS core. An absolute energy 
calibration of the EAS atmospheric Cherenkov light flux has been obtained. 
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1. Experiment QUEST and Simulations 

The QUEST experiment was developed to combine atmospheric Cherenkov light data 
to the charged particle EAS-TOP measurements (Gran Sasso, Italy, 2000 m a.s.l.) 1 The 
Cherenkov light detector was based upon five wide angle QUASAR-370 (37 cm diameter) 
semispheric photomultipliers installed on five telescopes (average pointing at the direction 
= 34°, ip = 167°, full field of view for 6 < 40°, used in analysis, is 0.41 sr). 

A new fitting function has been used to derive two main parameters of EAS Cherenkov 
light LDF for every recorded event: the light flux at core distance of 175 m Q175 and 
the LDF steepness, defined as the ratio of the fluxes at 100 and 200 m from the axis: 1 
P = Q(100)/Q(200). 

The energy measurement methods are based on parameterizations of Cherenkov light 
lateral distributions obtained from simulations performed through CORSIKA/QGSJET 2 "* 
for primary protons and iron nuclei with total energy 1,2,4 and 8 PeV, and zenith angles 9 
from 0° to 39°. 

The experiment has been additionally simulated by using such parameterizations and 
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Fig. 1. CORSIKA: Correlation between EAS 
Cherenkov LDF steepness and N e /Eg . 



Fig. 2. EXPERIMENT: Comparison of energy, ob- 
tained with two different methods for Echeb > 
2 X 10 15 eV. 



extracting the primary energy, the depth of EAS maximum X max , shower axis direc- 
tion and coordinates. The primary composition is assumed 50% p and 50% Fe. All other 
parameters have been calculated as functions of the main (i.e. P = F± (X ma ^,9, A), 
N e = F 2 (Eq,P,A) and so on). All the known apparatus errors and fluctuations have 
been taken into account. 



2. 1st Method: Size and Cherenkov Light LDF steepness P 

To derive the size comparable to the experimental one, we have taken into account both 
electrons and muons and used the experimental procedure of size reconstruction with 
NKG fitting function. Figure 1 shows the correlation between the parameter P (P = 
Q(100)/Q(200)) and the ratio of the size to primary energy (N e /Eo). Using this relation 
we can get the primary energy from measurement of N e and P: 

Esize [eV] = 1.59 x 10 11 iV e /10°- 33P . 

This method is quite composition independent, since it gives only 6% of difference in 
energy estimation for proton and iron. Its main advantage relies in the well developed 
technique of scintillator response calibration based on the measurement of the single par- 
ticle response. 5 From the simulation we derive for the reconstructed energy relative errors 
(mostly due to apparatus error of parameter P reconstruction) of about 30%. 

Similar method of energy reconstruction, but for LDF steepness, estimated at smaller 
distances from the core (20 — 100 m), was suggested in Ref. 6. 
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3. 2nd Method: Cherenkov Light Flux at 175 m Core Distance Q175. 

From the simulation we obtain the following relation between Q 175 [photon • cm" 2 • eV" 1 ] 
and primary energy Echer [eV]: 

Echer = 4 x 10 14 Q° 17 f. 

This method leads to larger differences in energy estimations for primary proton and iron 
(about 19%), than the first one. But the error of reconstructed energy (estimated with the 
above mentioned simulation of the experiment) is about 15% even for the assumed "50% p 
+ 50% Fe" composition, i.e. two times better than the first one. 

The main problem of this method is in the absolute calibration of Cherenkov light de- 
tectors, its error being estimated between 18% to 30% for different experiments. To get 
better accuracy we used the most probable ratio Esize / Echer as the coefficient for 
absolute calibration of Cherenkov detectors response. The coefficient is estimated for the 
natural mass composition. After such recalculation, the final comparison of two methods 
in the experimental data is shown in Fig. 2. The standard deviation of the experimen- 
tal distribution is very close to the one obtained in MC simulation of experiment for the 
SIZE/CLDF method (see above). 

4. A Reference Integral Primary Cosmic Rays Intensity 

600 events have been recorded during 140 h of data taking with energy (derived with the 
"re-calibrated Cherenkov light flux method") larger than 3 x 10 15 eV. The corresponding 
integral intensity is: 

/(£ >3x 10 15 eV) = (2.3 ± 0.1) x 10~ 7 rrr 2 • s" 1 • ster" 1 . 

The statistical uncertainty of the measurement is about 5%. 

This estimation can be used as a reference point for other experiments having no precise 
absolute calibration, such as Cherenkov light experiment TUNKA. 7 
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